Gravitational Wave Emission from a Bounded Source: the Nonlinear Regime 
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. We study the dynamics of a bounded gravitational collapsing configuration emitting gravitational 

' waves, where the exterior spacetime is described by Robinson- Trautman geometries. The full non- 

^S) ' linear regime is examined by using the Galerkin method that allows us to reduce the equations 

governing the dynamics to a finite-dimensional dynamical system, after a proper truncation proce- 
■ dure. Amongst the obtained results of the nonlinear evolution, one of the most impressive is the 

fact that the distribution of the mass fraction extracted by gravitational wave emission satisfies the 
distribution law of nonextensive statistics and this result is independent of the initial configurations 
considered. 
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Robinson- Trautman (RT) metrics are the simplest known solutions of vacuum Einstein's equations which 
may be interpreted as representing an isolated gravitational radiating systemP] . The metric can be expressed 
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■ where r is an afHne parameter defined along the shear-free null geodesies determined by the vector field 

qh, d/dr, dot means derivative with respect to u] \{u, 9) = — + "F^ ~ ^ cot 9 is the Gaussian curvature 

of the surfaces {u = const, r = const.). The geometry is nonstationary and axially symmetric, admitting 
the obvious Killing vector d jdip. As a consequence of Einstein equations, we have the following evolution 
equation 



In the above the subscript Q denotes derivatives with respect to 9. Eq. |5J is denoted the RT equation, 
governs the dynamics of the gravitational field and will be the basis of our analysis of the gravitational wave 
emission processes in RT spacetimes. Formally speaking, it allows to evolve initial data K(u, 9) prescribed 
on a given null surface u — uq (except in the case mg = 0). A particular and important solution of the field 
equations is the Schwarzschild metric obtained when K{u,9) = Kq = const, where X{u,9) = Kq^ and the 
corresponding Schwarzschild mass is Mschw = "t-o^o- This expression will be of particular importance in 
our characterization of the mass function of the configuration. Concerning this point, we note that in general 
the function M{u,9) = B{u)K^(u,9) is invariant under the coordinate transformation that reduces B{u) to 
a constant —2mo; also this transformation induces that the type-D 0{l/r^) curvature scalar associated with 
the mass aspect of the spacetime has a {u,9) dependence given exactly by BK^. 

Foster and Newman,^] exhibited the linearized solution of the field equations. We generalize their approach 
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in a rather distinct way in order to take into account the nonhnearities of the field equations. For this pro- 
posal, we shall consider the Galerkin projection method 4] in which we adopted the following decomposition 



K\u, x) = Al e«("'-) = Al exp ^ h,,{u)Pu{x) , (3) 

where we have introduced a new variable x — cos 9 such that —l<x< 1; iVis the order of the truncation, 
bk{u) are the modal coefficients, and we have adopted the Legendre polynomials Pk{x) as the basis functions 
of the projective abstract space. The internal product defined in this abstract space is {Pj{x), Pk{x)) = 
J^-^ Pj{x)Pk{x) dx — 2k+i ■ "^^"^ evolution equations for the modal coefficients bk{u) are derived in following 
way: we substitute (|3J| into the equation for A(u,x), both into |2Jl, and finally projecting the resulting 
equation into each basis function Pk(x). Then, 



bn{u) = 



(2n + l) / 



12moA^ 



[{l-X^)\']',Pn{x) 



(4) 



with n = 0, 1, 2.., N . Therefore, the dynamics of RT spacetimes is reduced to this system of (A'^+l) equations. 

The initial conditions bk{uo) are determined by the initial data K(uq, x) = k{x) since bj{uo) ~ —^-jp^-p^y^^ 
j = 0,..,N. We integrate numerically the set of equations for distinct families k{x) and the fol- 
lowing basic features of the nonlinear dynamics of RT spacetimes emerged, (i) All modal coefficients 
tend to zero asymptotically; the exception is bo{u) that approaches to a constant value which we denote 
by 6o(oo). This asymptotic configuration actually corresponds to the Schwarzschild solution with mass 
Moo = "IqAq exp(36o(oo)/2). Indeed, such a feature feature is in agreement with previous analytical stud- 
ies on RT spacetimes 0- (ii) For small values of the initial mass of the distribution, Minit, more rapidly 
the Schwarzschild solution is settled down; in addition a large fraction of Minit is lost in the process due 
gravitational wave extraction. On the opposite side, if Minn is large, then the evolution of RT spacetime 
towards the asymptotic Schwarzschild final state is very slow. Correspondingly, the fraction of Mi„it lost is 
very small, (iii) The invariant characterization of the presence of gravitational radiation can be done using 
the Peeling theorem. From this the complete pattern gravitational wave emission is obtained^- (iv) Finally, 
the last important feature is the relation between the fraction of mass extracted A = versus the 

initial mass Minit- For any family k{x) the mentioned relation satisfies the distribution provided by the 
nonextensive statistics 0. 




FIG. 1: Log-linear plot of A versus Minit- The solid line corresponds the nonextensive distribution A = Co (yo ~ 
y)" [1 + {q — 1)A2 (yo — y being associated to the initial mass, whereas the points were generated after 

integrating the dynamical system 0. The best fit is provided since yo = 1, Co — 4.462, a ~ 1.992 and q ~ 1.725. 
The initial data are ellipses with distinct eccentricities. 
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